Trapped p-wave superfluids: a local density approach 
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The local density approximation is used to study the ground state superfluid properties of har- 
monically trapped p-wave Fermi gases as a function of fermion-fermion attraction strength. While 
the density distribution is bimodal on the weakly attracting BCS side, it becomes unimodal with 
increasing attraction and saturates towards the BEC side. This non-monotonic evolution is related 
to the topological gapless to gapped phase transition, and may be observed via radio-frequency 
spectroscopy since quasi-particle transfer current requires a finite threshold only on the BEC side. 
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Recent experiments measuring momentum distribu- 
tion, collective modes, order parameter, quantized vor- 
tices, etc. have provided strong evidence for observa- 
tion of a superfluid phase in two-component cold atomic 
mixtures, interacting with short-range attractive s-wave 
interactions [E S U 11, IE @1 ■ These experiments have 
also shown evidence that the ground state of these s- 
wave mixtures evolves smoothly from a paired Bardeen- 
Cooper-Schrieffer (BCS) superfluid to a molecular Bose- 
Einstein condensate (BEC) as the attractive interac- 
tion varies from weak to strong values, marking the 
first demonstration of theoretically predicted BCS-BEC 
crossover 0,11,11]. 

On the other hand, there has also been substantial 
experimental progress studying p-wave interactions to 
observe triplet superfluidity il0|, [H [Hill, El ■ How- 
ever, controlling the p-wave interactions is proving much 
more difRcult due to the narrow nature of p-wave Fes- 
hbach resonances as well as to the more dramatic two- 
and also three-body losses [13, HH, E [ll [3- These 
experiments still motivated considerable theoretical in- 
terest predictin g q uantum and topological phase tran- 
sitions [H [H, IH m, E HI. More recently, p-wave 
molecules have been produced and their two-body prop- 
erties have been studied [2l|, opening the possibility of 
studying many-body properties of p-wave superfluids in 
the near future. 

In this manuscript, unlike the previous works on ho- 
mogenous systems [H, [H, [O, H, [H, HI], we study 
the ground state superfluid properties of harmonically 
trapped p-wave Fermi gases as a function of fermion- 
fermion attraction strength. Our main results are as 
follows. While we find that the density distribution is 
bimodal on the weakly attracting BCS side where the lo- 
cal chemical potentials are positive everywhere inside the 
trap, it becomes unimodal with increasing attraction and 
saturates towards the BEC side where the local chem- 
ical potentials become negative. This non-monotonic 
evolution is related to the topological gapless to gapped 
phase transition occurring in p-wave superfluids, and is in 
sharp contrast with the s-wave case where the superfluid 
phase is always gapped leading to a smooth crossover. 
Lastly, we propose that the phase transition found in the 



p-wave case may be observed via radio-frequency spec- 
troscopy since quasi-particle transfer current requires a 
finite threshold only on the BEC side, which is in sharp 
contrast with the crossover physics found in the s-wave 
case where a finite threshold is required throughout BCS- 
BEC evolution. 

Local density (LD) approach: To obtain these results, 
we consider a harmonic trapping potential, separate the 
relative motion from the center-of-mass one, and use LD 
approximation to describe the latter. In this approxi- 
mation, the system is treated as locally homogenous at 
every position inside the trap, and it is valid as long as 
the number of fermions is large which is typically satis- 
fied in cold atomic systems [22| ■ In order to describe the 
pairing correlations occurring in the relative coordinates, 
we use the BCS mean field (MF) formalism and neglect 
fiuctuations. The MF description is qualitatively valid 
throughout BCS-BEC evolution only at the low tempera- 
tures considered here ^,^8, 9], and it has been extensively 
applied to cold atomic systems describing qualitatively 
the experimental observations. 

Therefore, we start with the following local MF Hamil- 
tonian (in units of 7i = /cg 1) 
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where £ = {£ — 1) corresponds to s-wave (p-wave) 
systems, aj^ ^ creates a pseudo-spin a fermion with mo- 
mentum k, ^^(r, k) = e(k) — fj.i(r) is the dispersion with 
e(k) = k'^/{2m) and ^l|-}r) = - V{r). While the global 
chemical potential /i£ fixes the total number of fermions, 
the local chemical potential ^i{r) includes the trapping 



potential V{r) — 



V2. In Eq. A^(r,k) 



A^(r)r£(k) is the local MF order parameter where 
^eif) = 5 Ff (k)(a_k,|ak4) describes the spatial de- 
pendence, such that g > is the strengh of the attractive 
fermion-fermion interactions and (...) implies a thermal 
average. F£(k) determines the symmetry of the order pa- 
rameter given by F^(k) = W^(fc) Af,„i>f,m(k) where 
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W£{k) — k^knjlk^ + /co)(^+^^/^ describes the momentum 
dependence [Tl,[3- Here, ko ~ R^^ sets the momentum 
scale where Rq is the interaction range in real space. In 
this manuscript, we assume A^(r) to be real, and consider 
only the Ai.o = 1 and Ai_±i = symmetry. However, 
we emphasize that our discussion applies qualitatively to 
other p-wave symmetries as well. 

The local MF Hamiltonian can be solved by using stan- 
dard techniques [H, [l^, leading to a set of nonlinear 
equations for Ae(r) and fJ.i{r). These equations are 
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where ai is the experimentally relevant scattering param- 
eter which regularizes g, Et{r,'k) = ^^^^ (r, k) -I- (r, k) 
is the local quasi-particle energy, T is the temperature 
and V is the volume. In Eq. (jS]), nt{r) is the local density 
of fcrmions, and the total number of fcrmions N is fixed 
by = J drniir). Notice that ag has units of length 
(volume) in the s-wave (p-wave) case, and that our self- 
consistent solutions also describe the single pseudo-spin 
p-wave systems (except for the a summations here and 
throughout), which are presented next. 

BCS-BEC evolution in homogenous systems: To un- 
derstand the ground state properties of harmonically 
trapped p-wave superfluids within the LD approach, it 
is very useful to analyze first the homogenous s- and 
p-wave systems where V{r) = 0. Thus, next we dis- 
cuss the s-wave case where As(k) = A^ Wo(fc)^o,o(k) 
with yo,o(k) = 1/V47r, and compare these results with 
the p-wave case where Ap(k) = ApWi(A:)Yi^o(k) with 
Yi_o(k) = ^3/(47r) cos(6'k). In the numerical calcula- 
tions, while we mainly consider /cq — 100A:f to describe 
realistically the short-ranged atomic interactions, some 
of the fco = IQkp results are also shown for comparison. 

In Fig. [TJa), we show Ag and at zero tempera- 
ture (T = 0) for the s-wave case where the BCS-BEC 
evolution range in \/{kFas) is of order 1. Notice that. 
As grows continuously without saturation with increas- 
ing attraction, while /is decreases continuously from the 
Fermi energy ep = kp/{2m) on the BCS side to the half 
of the binding energy th,sl2 = —l/{2ma'^) on the BEC 
side [oj. Here, kp is the Fermi momentum which fixes 
the total density n = kp/{6TT'^) of fcrmions. Thus, 
we conclude that the evolution of Aj and jJi s 0jS 0, func- 
tion of l/{kpas) is analytic throughout, and BCS-BEC 
evolution is a smooth crossover 7, 8, 9]. 

In Fig. [TJb), we show Ap and /ip at T = for the 
p-wave case where the BCS-BEC evolution range in 
l/ikokpttp) is of order 1. Notice that, Ap is exponentially 
small but still finite in the BCS limit when fip ep and it 
grows rapidly with increasing attraction but almost sat- 
urates for large l/(fcofc|,ap), while /ip decreases continu- 
ously from ep on the BCS side to eb,p/2 = — l/(mfcoap) 
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FIG. 1: We show (in units of ep) the chemical potential fie 
(left y-axis) and the amplitude of the order parameter Ae 
(right y-a^xis) for (a) s-wave systems as a function of l/ikpas), 
and (b) p-wave systems as a function of l/{kok%ap). Here, 
solid (dotted) lines corresponds to fco — lOOkp (ko — 10^^)- 



on the BEC side. However, both Ap and fip are non- 
analytic exactly when /ip = at l/(fc|,ap) w 0.45, which 
occurs on the BEC side of unitarity. We note that the 
non-analyticity of /ip is barely seen in Fig. (Hb), and it 
is more expilicit in derivatives of /Xp. Thus, in the p- 
wave case, BCS-BEC evolution is not a crossover, but a 
quantum phase transition occurs [H, [13, ■ 

This phase transition can be understood as follows. 
The quasi-particle excitation spectrum EiCk) is gapless 
when the conditions Af(k) = and £,e{k.) = are both 
satisfied for some k-space regions. While the second con- 
dition is satisfied for both s- and p-wave symmetries on 
the BCS side where /i^ > 0, the first condition is only 
satisfied by the p-wave order parameter. Therefore, un- 
like the s-wave case, i?p(k) is gapless on the BCS side 
(//p > 0) but it is gapped on the BEC side (up < 0), lead- 
ing to the phase transition discussed above |3. [isl. UM. [23| . 

Having discussed the ground state of homogenous sys- 
tems, next we analyze the trapped case. 

BCS-BEC evolution in trapped systems: For this 
purpose, similar to the analysis of homogenous sys- 
tems, first we discuss the s-wave case where As(r, k) = 
As{r)Wo{k)Yofi{'k), and compare these results with the 
p-wave case where Ap(r, k) = Ap{r)Wi{k)Yifi{'k). In 
the numerical calculations, we again choose ka = lOOkp 
where kp — mujorp is the global Fermi momentum. 
Here, rp is the Thomas-Fermi radius determined by 
V{rp) = ep = kf/(2m), and fixes the total number of 
fcrmions to = J2a- k%r'p/AS,. 

Within the LD approximation, the density distribution 
of trapped non-interacting (5 = or ag — )■ 0^) gas is 
T^iif) = X^cr ^f(^)/(6i'^); whcrc kp{r) is the local Fermi 
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FIG. 2: (Color online) We show [in units of fc|,/(27r)^] the 
density distribution n^(r) for (a) s-, and (b) p-wave systems 
as a function of the trap radius r (in units of rp)- Here, we 
set ko = lOOfcF. 



best on the BEC side where strongly attracting fermion 
pairs form weakly repulsive local molecules which can be 
well described by the Bogoliubov theory 9, 19]. On this 
side, the size of the s-wave molecules decreases to arbi- 
trarilly small values as S.b,s as > when fcoa^ ^ 1, 
leading to arbitrarilly weak molecule-molecule repulsion 
Ubb.s = ZTrasB.s/Ti where qbb.s = 2as within the 
Born approximation Q. However, the size of the p- 
wave molecules saturates to small but finite values as 
^B,p ~ when fegOp ^ 1, leading also to a weak 

but finite molecular repulsion U bb,p = 27raBB^p/m where 
dBB.p — Q/ko [l^. Since our LD approximation recov- 
ers the Thomas-Fermi approximation for the resultant 
molecules in the BEC limit, Tipir) saturates rapidly for 
l/{kokpap) > due to the presence of weak but finite 
molecular repulsion. 




(a) BCS (b) Unitarity (c) BEC 



momentum determined by fii = kp{r) / {2m) + V{r) with 
lie = ep- Therefore, both kpir) and ne{r) are highest at 
the center of the trap as can be also seen in Figs. [D^a) 
and[2l[b) when l/(kpas) — — oo and l/^k^kpap) — —oo, 
respectively. 

In the presence of weak attraction, while fie deviates 
from ep, the density distribution is still well described 
by the non-interacting expression. For fixed N, ne{r) 
is expected to squeeze and increase towards the center 
of the trap since /i^ decreases with increasing attrac- 
tion. The squeezing effect of the weak attractive inter- 
actions can be seen in Fig. for the s-wave and in 
Fig. [2tb) for the p-wave systems when l/ikpas) = —0.5 
and l/{kokpap) — —0.4, respectively. 

However, weakly attracting p-wave superfluids have bi- 
modal density distribution as shown in Fig. [2Jb), which 
is in sharp contrast with the unimodal s-wave distribu- 
tion of Fig. mja). This difference can be understood from 
the homogenous results shown in Fig. [Tfa) and [2b) as 
follows. First, since fci?(r) decreases away from the center 
of the trap, the local s- and p-wave scattering parame- 
ters l/[kp{r)as] and l/[kQkp{r)ap], respectively, increase 
as a function of r if and ap are fixed. Second, notice 
in Fig.[TJb) that Ap increases rapidly from exponentially 
small to larger values as a function of 1 / {k^k'j^ap) , un- 
like As which increases smoothly as shown in Fig. HJa). 
These two observations combined shows that the almost 
non-interacting ripir) distribution towards the tail is due 
to finite but exponentially small Ap(r). 

With increasing attraction towards unitarity, while the 
unimodal ns(r) distribution smoothly squeezes further as 
shown in Fig. [2ja) for l/{kpas) = and l/{kpas) = 0.5, 
the bimodal np{r) distribution becomes unimodal and 
saturates as shown in Fig. [2)Jb) for l/{kokpap) = and 
l/{kQkpap) — 0.4. This difference can be understood at 



FIG. 3: Schematic diagrams showing (a) a fully gapless su- 
perfluid on the BCS side when ^p(0 < r ^ rp) > 0, (b) 
a partially gapped superfluid around unitarity when fip{Q < 
r < r*) > but fip{r > r,) < 0, and (c) a fully gapped 
superfluid on the BEC side when fip{r > 0) < 0. 

This non-monotonic evolution of np{r) is also related 
to the topological phase transition discussed above for 
the homogenous systems. In the trapped case, the lo- 
cal quasi-particle excitation spectrum i?£(r, k) at posi- 
tion r is gapless when the conditions A^(r, k) = and 
^^(r, k) = are both satisfied for some k-space regions. 
While these conditions are both satisfied everywhere in- 
side the trap leading to a fully gapless superfluid on the 
BCS side, they are only satisfied around the center of 
the trap close to unitarity leading to a partially gapped 
superfiuid. Further increasing the attraction towards the 
BEC limit, the second condition is not satisfied, and the 
entire superfiuid becomes fully gapped. These phases are 
schematically shown in Fig. [n^a), [3Jb) and[3]Jc), respec- 
tively, and next we discuss their experimental detection. 

Radio-frequency (RF) spectroscopy: The gapless to 
gapped phase transition discussed above may be observed 
for the first time in cold atomic systems via for instance 
RF spectroscopy, where atoms are transferred from one 
hyperfine state to another generating a quasi-particle cur- 
rent [22, 24, 25, 26]. This is analogous to electrons tun- 
neling from a superconducting to normal metal, and it 
has been used in atomic systems to observe pairing cor- 
relations in unpolarized [2^ as well as polarized [1^ mix- 
tures. 

The local quasi-particle transfer current, within the LD 
approximation, is given by [13, 113, 113 

h{r, co)^t%Y, MK Ur, k) - oo]F[Ur, k) - a;], (4) 
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FIG. 4: (Color online) We show (in units of pFt% /2) the quasi- 
particle transfer current Ie{io) for homogenous (a) s-, and (b) 
p-wave systems as a function of the effective detuning (jj (in 
units of ep)- Here, we set fco = lOOfcF- 

where tp is the transfer amplitude, Ai(\i,x) is the spec- 
tral function corresponding to the superfluid state, and 
F{x) — l/[exp(a:/T)-|-l] is the Fermi function. Here, uj — 
ujl — ojh is the effective detuning where lol and loh are 
RF laser frequency and hyperfine splitting, respectively. 
We evaluate Eq. (|4|) with the standard BCS spectral func- 
tions ^^(k, e) = 27r{u2(r, k)^[e - Ei{r, k)] + v]{r, k)S[e + 
Eeir,^)]}, where u'j{r,k) = 0.5[1 + ^i{r,k)/ Ee{r,k)] and 
Vg{r, k) = 0.5 [1 — (r, k)/i?^(r, k)] are coherence factors, 
and S{x) is the delta function. 

At T = 0, the s-wave current /^(r, w) can 
be evaluated analytically leading to /s(r, w) = 
npFtl [As {r)/{4TTu;)]^ y'C{r,uj)e{ij)9[C{r, cj)] , where 
Pp — mVkF/{2TT^) is the density of states. 



C{r,uj) = [lu^ ~ A2(r)/(47r)]/(2tj) + ^,(r), and e{x) is 
the theta function. Therefore, Is{r,u!) flows when the 
threshold LOthAr) > -M.(^) + v^MsCO + A2(r)/(47r) > 
is reached, which reduces to ujth,s{r) > A^(r)/[87r^s(»^)] 
in the BCS and ujth,s{r) > 2\p.s{r)\ in the BEC limit. 
Therefore, 'Wth,si'r) 7^ everywhere inside the trap 
throughout BC'S-BEC evolution [H H, HI, . These 
finite detuning thresholds can be also seen in Fig. [H^a), 
where we show /s(a^) for homogenous systems. 

The p-wave current /p(r, w) is difficult to evaluate ana- 
lytically. However, unlike the s-wave case, we expect that 
Wth,p{r) 7^ ^ everywhere inside the trap only on the BEC 
side beyond unitarity where Ep{r,k) is gapped, and also 
that wth.p{r) = everywhere inside the trap on the BCS 
side where Ep{r, k) is gapless. The absence (presence) 
of finite thresholds in gapless (gapped) superfluids can 
be seen in our homogenous results shown in Fig. IH^b). 
Notice that, while ujth,p = for l/{kokpap) — —00, —0.4 
and 0, a finite threshold is required for l/{kokpap) — 0.4. 
Based on these results, we hope that spatially resolved 
RF spectroscopy measurements (similar to [26|) may be 
used to identify all three phases proposed in Fig. [31 

Conclusions: To summarize, we showed that while the 
density distribution of p-wave systems is bimodal on the 
weakly attracting BCS side, it saturates and becomes 
unimodal with increasing attraction towards the BEC 
side. We discussed that this non-monotonic evolution is 
related to the topological gapless to gapped phase tran- 
sition occurring in p-wave superfiuids, and is in sharp 
contrast with the s-wave case where the superfluid phase 
is always gapped leading to a smooth crossover. Lastly, 
we proposed that this phase transition may be observed 
via RF spectroscopy since quasi-particle transfer current 
requires a finite threshold only on the BEC side, which 
is in sharp contrast with the s-wave case where a finite 
threshold is required throughout BCS-BEC evolution. 
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